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PREFACE 

This  report  was  prepared  by  Dr.  Motoi  Kianai,  Research  Physicist,  Research 
Division,  Cold  Regions  Research  and  Engineering  Laboratory  (CRREL).  U.S.  Army 
Terrestrial  Sciences  Center  (USA  TSC).  Computer  work  was  done  'ey  SP5  Jack  D. 
Russell  of  CRREL. 

The  report  covers  the  optical  and  physical  studies  of  ire  fog  undertaken  by 
CRREL  as  a  FY68  ILIR  project  with  partial  support  from  the  U.S.  Army  Arctic  Test 
Center.  This  report  was  published  under  DA  Project  lT06r*01A91A.  lo-Honst  Lab¬ 
oratory  independent  Research. 

USA  TSC  is  a  research  activity  of  the  Army  Mattel  Command. 
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ABSTRACT 


Ice-fog  crystals  consisting  of  many  spherical  particles,  and  some  hexagonal 
plates  and  columns,  were  observed  at  ambient  temperatures  of  about  -IOC  in  the 
Fairbanks.  Alaska,  area  during  mid-winter.  The  concentrations  and  the  size 
distributions  of  the  ice-fog  crystals  were  measured.  The  attenuation  and  back- 
scattering  of  infrared  radiation  by  ice-fog  crystals  were  computed  for  optical 
wavelengths  of  2.2ft,  2.7^,  4.fy,  5.75^,  9.7^1  and  10.9^  using  the  Mie  theory.  The 
minimum  attenuation  coefficients  and  backscattering  functions  of  ice  fog  were 
found  to  be  at  9.7pi  wavelength  in  the  observed  wavelengths.  Optical  attenuation 
coefficients  and  volume  backscattering  functions  of  water  fogs  were  also  computed 
using  the  Mie  theory.  The  minimum  attenuation  coefficients  and  backscattering 
functions  of  water  fog  were  found  to  be  at  10.9^  wavelength  in  the  region  of  2.2^, 
2.7 ft,  4.5 fi,  5.75ft,  9.7 y.  and  10.9^.  Both  the  attenuation  coefficients  and  back- 
scattering  functions  of  ice  fog  are  within  the  same  order  of  magnitude  as  water 
fog  for  equivalent  fog  concentrations  and  wavelengths. 


THE  ATTENUATION  AND  BACKSCATTERING  OP  INFRARED 
RADIATION  BY  ICE  FOG  AND  WATER  FOG 

by 

Motoi  Kumai  and  Jack  D.  Russell 


INTRODUCTION 

By  reducing  visibility,  ice  fog  hampers  the  normal  flow  of  traffic  during  severe  cold  weather 
conditions  in  subarctic  countries.  In  addition,  ice  fog  can  cause  extinction  of  the  infrared  beam 
in  an  infrared  guidance  system.  Because  of  this,  the  emphasis  of  CRREL's  ice  fog  research  is 
being  placed  on  the  effects  of  ice  fog  on  infrared  transmission.  This  report  presents  results  of  a 
study  of  attenuation  and  backscaltering  of  infrared  radiation  by  fog  in  Alaska. 


PHYSICAL  PROPERTIES  OF  ICE  FOG 

The  physical  properties  of  ice  fog  in  the  Fairbanks,  Alaska,  area  have  been  studied  by  Thuman 
and  Robinson  (1954).  The  formation  of  ice  fog  and  its  nuclei  have  been  studied  by  Kumai  (1964) 
and  Kumai  and  O'Brien  (1965). 

The  optical  properties  of  fog  depend  on  the  number  concentration  and  size  distribution  of  the 
particles,  which  can  vary  significantly  during  different  meteorological  conditions.  Establishing 
a.  representative  fog  model  is  not  an  easy  task.  Determination  of  the  collection  efficiency  for  fog 
droplets  with  radii  approaching  the  wavelength  of  near-infrared  radiation  is  one  of  the  problems  in 
fog  sampling. 

Ice-fog  crystals  and  ice  crystals  are  initial  stages  in  the  formation  of  snow  crystals.  Be¬ 
cause  their  sizes  and  optical  properties  are  quite  different,  it  is  convenient  to  consider  them 
separately.  Ice  crystals  have  well-formed  hexagonal  plates  and  columns.  They  are  often  called 
“diamond  dust"  because  of  the  way  they  twinkle  under  a  lew  sun  or  under  nighttime  illumination. 

Ice  crystals  form  at  about  -25C  and  are  20^  to  309^  in  diameter.  They  have  much  lower  concentra¬ 
tions  and  do  not  affect  visibility  as  adversely  as  ice  fog  does.  Ice-fog  crystals  consist  of  many 
spherical  particles,  equant  crystals  with  rudimentary  faces,  and  some  hexagonal  plates  and  columns 
of  2-ji  to  30  (i  diameter  formed  at  ambient  temperatures  of  about  —10C  by  spontaneous  nueleation 
(Fig.  1). 

Ice-fog  distributions  at  air  temperatures  of  -39C  and  -4 1C  in  the  Fairbanks,  Alaska,  area  are 
shown  in  Figure  2  and  Figure  3.  The  number  of  particles  and  the  mass  of  fog  per  unit  volume  are 
shown  as  a  function  of  crystal  diameter.  Table  1  summarizes  the  important  physical  properties  of 
each  fog  spectrum. 
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Figure  3.  Size  and  miss  distribution  ol  ice-log  crystals  formed  at 
-41C  ambient  'emperature. 


Table  I.  Physical  properties  of  ice  fog  at  Fairbanks,  Alaska. 


Observations 

N 

rmode 

rmin 

rmax 

Ar 

Air  temp 

L.W.C. 

(ao. /cm1) 

<P> 

<Fl 

(P> 

<P> 

(°C) 

No.  1  (Fig.  2) 

140 

3.0 

1,5 

12,0 

n 

-39 

0.08 

No.  2  ( Ftg.  3) 

90 

1.5 

*1.5 

12.0 

HI 

-4 1 

0.02 

N  =  total  concentration 

rmode  ~  mode  radius  =  radius  corresponding  to  the  maximum  number  of  ice-fog  crystals 
rmin  =  m'nitn'im  t&dius 

r_„.  =  maximum  radius 

max 

Ar  =  radius  interval  containing  r.(r)  crystals,  where  N  =  £  n( r) 

r 


L.V.C.  =  liquid  water  content. 
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MIE  SCATTERING  COMPUTATIONS 

The  Mie  theory  is  ?.a  exact  solution  of  Maxwell’s  electromagnetic  theory  for  spheres;  it 
describes  the  field  inside  and  outside  a  spheie  for  any  scattering  angle,  size,  and  index  of  refrac¬ 
tion. 

For  backscattering,  the  scattering  angle  equals  -  radians.  The  particle  size  parameter  is 
x  -  k' a  where  k‘  =  2jt/A.  a  is  the  radius  of  the  spheres,  and  A  is  the  wavelength  of  incident  radia¬ 
tion. 

When  dealing  with  actual  particles,  scattering  due  to  the  sphere’s  dielectric  properties  alone 
is  not  sufficient  to  define  fully  the  intensity  of  the  backscattered  radiation;  the  contribution  due 
to  absorption  within  the  sphere  must  also  be  included,  This  requires  a  complex  index  of  refrac¬ 
tion  of  the  form  m  =  n  -  iJc,  where  m  - (<  -  4 nio/tj)'1,  f  being  the  dielectric  constant,  o  the 
conductivity,  and  u  the  angular  frequency  of  the  electromagnetic  wave. 

The  computational  method  as  outlined  in  Carrier  et  al.  ( 1967)  appears  quite  simple  a:  first 
glance.  In  the  present  study  we  me  interested  in  the  attenuation  coefficient  due  to  scattering,  b, 
and  the  volume  backscattering  function,  p(n).  The  quantity  b  is  the  coefficient  in  the  relationship 

/  =  /0  exp  (-M)  (1) 

defining  the  decrease  in  light  intensity  while  passing  through  a  medium  of  length  /. 

According  to  the  above-mentioned  paper, 

rmax 

b  =  2  n(r)8r/rr  Kgxl  (x,m)  (meter'1]  (2) 

rmin 

and 

A2  rsnax 

Bin)  =  — -  1  nfrlSri,  (n,  x,  m)  (3) 

4rr  r  . 
min 

where  nfr)  is  the  number  of  particles  per  unit  volume  per  Sr  radius  interval,  Kgxl  is  the  total 
extinction  cross  section  of  one  particle,  and  i,  is  the  Mie  intensity  function  of  a  scattered  compo¬ 
nent  whose  electric  vector  is  either  perpendicular  or  parallel  to  the  plane  of  observation.  When 
6  =  7,\  then  i,  =  i3. 

The  problem  becomes  more  difficult,  however,  when  we  examine  the  work  by  Deirmendjian  and 
Clasen  (1962),  which  outlines  a  computational  scheme  to  derive  values  for  Kext  and  i,.  The  method 
involves  the  use  of  ordinary  Bessel  functions,  and  circular  and  hyperbolic  functions  with  complex 
argument.  Although  the  functions  may  be  reduced  to  fairly  simple  recursion  formulas,  the  presence 
of  complex  quantities  in  the  denominator  of  most  of  the  terms  requires  long  and  tedious  algebraic 
manipulation  to  make  the  problem  programmable.  It  must  also  be  kept  in  mind  that  a  computer 
deals  with  pure  numbers  and  that  the  real  and  the  imaginary  parts  of  a  complex  number  must  be 
dealt  with  individually. 

The  greatest  consumer  of  computer  real-time  is  the  evaluation  of  the  recursion  formulas  which 
consist  of  very  slowly  converging  infinite  series.  It  is  not  unusual  to  find  that  300-350  terms  are 
required  to  solve  for  just  one  value  of  the  quantities  Kgxt  and  i,  of  which  there  are  as  many  values 
as  there  are  radius  intervals.  We  must  then  multiply  the  number  of  radius  intervals  by  the  number 
of  wavelengths  in  question  when  estimating  total  run-time. 
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In  Figures  £  and  3,  the  radius  interval,  Ar,  is  0.5//.  Particle  concentrations  for  each  fog 
radius  oirk)  were  punched  on  paper  tape  for  subsequent  entry  into  the  computer.  Computer  calcula¬ 
tions  of  the  attenuation  coefficient  and  backscattering  functions  from  eq  2  and  3  were  made  for 
the  fog  size  distributions  of  Figure  2  and  Figure  3,  using  the  optical  constants  of  ice  byKislovskii 
(19591  and  the  optical  constants  of  water  by  Centeno  (1941). 

The  recent  investigations  on  the  optical  properties  of  ice  were  done  by  Hornig  et  a l.  (1958), 
Ockraan  f  1958),  and  Zander  (1966).  However,  their  papers  do  not  present  the  values  of  the  refrac¬ 
tive  index  (n)  and  extinction  coefficient  (it)  of  ice.  Due  to  the  experimental  difficulties,  our 
present  knowledge  of  the  optical  constants  of  ice  in  the  infrared  region  is  inexact  compared  with 
our  knowledge  of  those  of  water.  The  optical  constants  of  ice  and  water  used  for  this  computer 
calculation  are  presented  in  Table  II  (Centeno,  1941;  Kislovskii.  1959). 


Table  II.  Optical  constants  of  ice  and  water. 


tee  _  _ Water 


Wavelength 

Extinction 

Refraction 

Extinction 

Refraction 

in  microns 

coelticient 

index 

coefficient 

index 

A 

k 

n 

k 

n 

2.2 

0.0005 

1.29 

0.0005 

1.293 

2.7 

0.03 

1.10 

0.0183 

1.216 

4.5 

0.01 

1.33 

0.0184 

1.341 

5.75 

0.01 

1.24 

0.0427 

1.273 

9.7 

0.08 

1.13 

0.0579 

1.230 

10.9 

0.20 

1.12 

0.0993 

1.150 

RESULTS 

Optical  attenuation  coefficients  and  volume  backscattering  functions  were  computed  for  the 
observed  Alaskan  ice  fogs  for  optical  wavelengths  of  2.2 g,  2.7//,  4.5//,  5,75//,  9.7ft  and  10.9ft  using 
the  Mie  theory.  The  calculations  were  made  for  ice-fog  crystal  concentrations  of  70,  140,  280  and 
420/cm3  with  size  distribution  shown  in  Figure  2;  and  ice-fcg  crystal  concentrations  of  90,  180, 
270  and  360/cm3  with  size  distribution  shown  in  Figure  3.  The  attenuation  coefficients  bfnr’)  and 
the  backscattering  functions  /3(ff)(m*,srI)  of  ice-fog  crystals  are  presented  in  Tables  III  and  IV. 
The  minimum  attenuation  coefficient  and  backscattering  function  are  found  to  be  at  9.7fx  wave¬ 
length  for  the  computed  wavelength  range  of  2.2//,  2.7//,  4.5/,  5,75//,  9.7//  and  10.9//.  The  optical 
wavelengths  of  2.2//,  9.7//  and  10.9//  are  those  corresponding  to  atmospheric  windows.  The  wave¬ 
lengths  of  2.7//  and  4.5//  fall  within  the  absorption  bands  of  both  water  vapor  and  carbon  dioxide, 
and  the  wavelength  of  5.75//  is  within  the  absorption  band  of  water  vapor. 

When  the  air  temperatures  were  above  -15C,  supercooled  fogs  wore  observed  in  the  Fairbanks, 
Alaska,  area  during  mid-winter.  Optical  attenuation  coefficients  and  volume  backscattering  func¬ 
tions  of  water  fogs  were  computed  for  optical  wavelengths  of  2.2//,  2.7//,  4.5//,  5.75//,  9.7 ft  and 
10.9//  using  the  Mie  theory.  The  calculations  were  done  for  fog-droplet  concentrations  of  70,  140, 
280  and  420/cm3  with  size  distributions  shown  in  Figure  2,  and  fog-droplet  conceni.ations  of  90, 
180,  270  and  360/cm3  with  size  distributions  shown  in  Figure  3.  The  attenuation  coefficients 
b(nr')  and  the  backscattering  functions  /3('rXm'Isr ')  of  water  fogs  are  presented  in  Tables  III  and 
IV.  Both  the  attenuation  coefficients  and  backscattering  functions  of  ice  fog  and  water  fog  are 
within  the  same  order  of  magnitude  for  the  same  fog  concentrations  and  wavelengths.  The 
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Table  m.  TV*  nix— iUn«  cooffldeaU  h( «'*). 


Sis*  Co*c*suttloa  Ltqsid  scalar 


itmrihnlo* 

ef  tog  dropiata 
ao./c*1 

CMtIK 

i/m  *  . 

M 

1.7 

4.5 

J7S 

17 

TO 

Ece-fo* 

Ft*,  e 

70 

0.039 

1. 163*  10“ 

1.2 16 10“ 

1.32  lx  10*1 

1.304x10“ 

5.136  UT 

6.795x10-* 

Ft*.  2 

140 

0.077 

2.306x10** 

2.424X10*1 

2.642x10“ 

2.606x10“ 

1.026x10“ 

1.356 1C“ 

Ft*.  2 

290 

0.15 

4.611x10“ 

4.848x10“ 

5.283x10“ 

5.21610“ 

2.052x10“ 

2.717x10“ 

Ft*.  2 

420 

0.23 

6.917x10“ 

7.272x10“ 

7.925x10“ 

7.824x10“ 

3.07610“ 

4.07610“ 

Ft*.  ! 

90 

0.023 

6.932x10“ 

6.694x10“ 

9.719x10“ 

7.936 10“ 

3.0*610“ 

4.30610“ 

<1*.  S 

180 

0.047 

1.796x10“ 

1.789x  10“ 

1.944x  10“ 

1.586 10*' 

6.092x10“ 

8.616(0“ 

Ft*.  3 

270 

0.070 

2.860X10“ 

2.606x10“ 

2.916<10“ 

2.362x10“ 

9.13610“ 

1.296  (O'* 

Ft*.  •* 

360 

0.093 

3. £73x10“ 

3.478x10“ 

im  io“ 

3.17610“ 

1.216 10“ 

1.723x10“ 

■alar-fo* 

FI*.  2 

70 

0.039 

1.153x10“ 

1.217x10“ 

1.314x10“ 

1.282x10** 

7.84610“ 

5.456x10“ 

Ft*.  2 

140 

0.077 

2.306x10“ 

2.434x  10“ 

2.626(1 0“ 

2.564x10'* 

1.568x10“ 

1.091x10“ 

FI*.  2 

290 

C.15 

4.611x10“ 

4.86610“ 

5.26610“ 

5.127x10“ 

3.137x10“ 

2.18610“ 

FI*.  2 

420 

0.23 

8.916(10“ 

7.803x10“ 

7.886x10“ 

7.89  U 10“ 

4.705x10“ 

3.273x10** 

FI*.  3 

90 

0.023 

6.919x10“ 

9.064x10“ 

9.875x10“ 

8.256x10“ 

4.504x10“ 

3.256x10“ 

FI*.  S 

no 

0.047 

1.784x10“ 

1.811x10“ 

1.935x10“ 

1.651x10“ 

9.006x10“ 

6.509x10“ 

FI*.  3 

270 

0.070 

2.873x10“ 

2.716xH>“ 

2.902x10“ 

2.4:7x10“ 

1.351x10** 

9.764x10“ 

FI*.  3 

360 

0.093 

3.666x10“ 

3.622x10“ 

1.870x10“ 

3.303x10“ 

1.602x10“ 

1.30610“ 

Table  IV.  Tbe  backer  marl  a*  (wUtu  flf.nKmr'wrt\ 


Slbe  CoacaatraUoa  Liquid  water 

dlstr&iiUoe  ot  log  dropiata  coaiaac  _ 

ao./cai* _ g/m' _ l.l 


Wtnhigth,  n 

l.l _ 4J _ SJS _ 9.7  10.9 


JeMo* 


Fl*.  2 

70 

0.039 

4.504x10“ 

9.18610“ 

1.251x10“ 

2.69610“ 

2/64xl0“ 

5.460x10'* 

Fl|.  2 

140 

0.077 

9.809x10“ 

1.817x10'' 

2.601x10“ 

5.389x10'* 

0.927x10“ 

1.090x10“ 

Fl*.  3 

280 

0.15 

1  92610“ 

3.654x10“ 

5.00610“ 

1.073x10“ 

1. 185x10“ 

2. 180x9)'* 

Fi*.  2 

420 

0.23 

2.88610“ 

5.481x10“ 

7.50610“ 

7.606 K  * 

1.773x10'* 

3.271x10“ 

Fl*.  8 

90 

0.023 

*.744x#“ 

5.794x10“ 

7. 5 16 10 '* 

1.936>  10“ 

2.486x10“ 

4.689x10“ 

Ft*.  3 

180 

0.047 

6. 489x10“ 

1.119x10“ 

1.50610“ 

8.37  *10“ 

4.970x10“ 

9.179x10“ 

Ft*.  3 

270 

0.070 

8.23610“ 

1.738x10“ 

2.25610“ 

5,PJ8xlO“ 

7.456x10“ 

1.377x10“ 

Fl*.  3 

880 

0.093 

1.098x10“ 

2,31Bx  10“ 

3.006x10“ 

7.745x10“ 

9.940x10“ 

1.886x10** 

Vater-fof 


Ft*.  2 

70 

0.039 

'  960x10“ 

*.699x10“ 

9.400-  10“ 

1.53610’* 

1.170xi0“ 

4. 10610“ 

Fl*.  2 

HO 

0.077 

9.900x10“ 

6.397x10“ 

1.880x10“ 

3.070x10“ 

HJSOxlO'* 

3.206x10“ 

Ft*.  2 

280 

0. 16 

1.996 10'* 

1.079x10“ 

8.730x10“ 

8. 18  lx  10“ 

4.679x10“ 

1.64  lx  10“ 

Ft*.  2 

420 

0.23 

2.988x10“ 

1.819x10“ 

6.640x10“ 

9.271x10“ 

7.018x10“ 

*46610“ 

Ft*.  3 

90 

0.023 

2.87610“ 

1.676 10'* 

5.990x10“ 

1.818x10“ 

8.98610“ 

8.466x85“ 

Ft*.  3 

180 

0.047 

6.145x10“ 

3.846x10" 

1.198x10“ 

2.63610“ 

1.78610'* 

6.961x10“ 

F%.  8 

270 

0.070 

8.618x10“ 

5.019x1'/* 

1.796x10“ 

3.948x10“ 

2.680x10“ 

L040xt0“ 

Ft*.  8 

860 

0.093 

1. 148x10“ 

6.89610'* 

2.89610“ 

5.285x10“ 

8.57610“ 

1.886x10“ 
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minimum  attenuation  coefficient  and  backscaltering  functions  of  water  fog  were  found  to  be  at 
10.9 fi  wavelength  for  the  observed  wavelengths  of  2.2ft.  2.7 ft,  4.5 ft,  5.75p,  9.7 ft,  and  10.  V 

These  calculations  desciibe  only  the  effect  of  water  droplets  and  ice  crystals  on  the  attenua¬ 
tion  and  scattering  of  radiation.  For  computations  of  total  attenuation  and  scattering,  the  effects 
of  the  atmospheric  absorption  bands,  especially  water  vapor  and  CO,,  would  also  have  to  be 
considered. 
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-Ice*-fog  crystals  consisting  of  many  spherical  paiticles,  and  some  hexagonal  plates 
and  columns,  were  observed  at  ambient  temperatures  of  about  -40C  in  the  Fxirbanki 
Alaska,  ,area  during  mid-winter.  The  concentrations  and  the  sise  distributions  of 
the  ice -fog  crystals  were  measured.  The  attenuation  and  backscattering  of  infrared 

radiation  by  iee-fog  crystals  were  computed  for  optical  wavelengths  of  2.  2n,  2.  7jjj _ 

4.  5p.  5.  75fi,  9.7p  and  10.  9p  using  the  Mie  theory.  The  minimum  attenuatioiTccef- 
ficients  and  backscattering  functions  of  ice  fog  were  found  to  be  at  9.  7p  wavelength 
in  the  observed  wavelengths.  Optical  attenuation  coefficients  and  volume  backseat  - 
tering  functions  of  water  fogs  were  also  computed  using  the  Mie  theory.  The  mini¬ 
mum  attenuation  coefficients  and  backscattering  functions  of  water  fog  were  found  to 
be  at  10.  9p  wavelength  in  the  region  of  2.  2p,  2.  7p,  4.  5p',"B.  75u,  9.  7p  and  10.  9p. 
Both  the  attenuation  coefficients  and  backscatteringi functions  of  ice  fog~are  within 
the  same  order  of  magnitude  as  water  fog  for  equivalent  fog  concentrations  and 
wavelengths.  ...i 
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